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Abstract 

We present a theory of tunneling spectroscopy of spin-selective Aharonov-Bohm oscillations 
in a lateral triple quantum dot molecule. The theory combines exact treatment of an isolated 
many-body system with the rate equation approach when the quantum dot molecule is weakly 
connected to the leads subject to arbitrary source-drain bias. The tunneling spectroscopy of the 
many-body complex is analyzed using the spectral functions of the system and applied to holes in 
a quantum dot molecule. Negative differential conductance is predicted and explained as a result 
of the redistribution of the spectral weight between transport channels. It is shown that different 
interference effects on singlet and triplet hole states in a magnetic field lead to spin-selective 
Aharonov-Bohm oscillations. 

PACS numbers: 73.21.La,73.23.Hk,73.63.Kv 
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I. INTRODUCTION 



Quantum dot (QD) systems are artificial quantum mechanical systems with great 
controllability of the electronic, spin, and transport properties.^ The occupancy and 
spins of electrons in single and coupled semiconductor QD systems can be controlled 
electrostatically.-'^'^'^'^'^ In lateral structures where QDs are defined by lateral gates on 
top of a heterojunction containing two-dimensional electron gas (2DEG), the confinement 
potentials and tunneling barriers can be easily tuned by applying voltages to the gates. Cou- 
pled QD systems, where the coherence of different quantum states can play an important 
part, are one of promising candidates for realizing quantum information and computation 
de vices. Coulomb interactions in such small systems become more important and QD 
networks offer great flexibility for studying many-body effects in a very controllable manner. 
Transport through a strongly interacting QD system leads to many interesting phenomena 
such as Coulomb blockade^ due to the strong Coulomb interaction and spin blockade in a 
double QD system with an electron localized in one of the dots.— High source-drain bias 
transport can be used as a spectroscopic tool for the interacting many-body system.— 

In triple quantum dot (TQD) systems, the electronic and spin properties depend on the 
topology of the system.— Interference effects are expected in triangular TQD molecules^^ 
and coherent transport was experimentally observed.-i^ We study a TQD system which is 
connected to two leads [Fig. [1] , with dot 1 connected to the left lead and dot 3 connected 
to the right lead. In linear response regime, for an electron to move from one lead to 
the other we need three different charge configurations {Ni, N2, N^), {Ni + 1,A^2,^3), and 
(A^^i, N2, + 1) in terms of the number of electrons in each dot to be on resonance. That 
is, E{Ni, N2, N3) + /io = E{Ni + 1, N2, N3) = E{Ni, N2, N3 + 1) where fio is the chemical 
potential of the leads. This allows an electron to move from the left lead to dot 1, to dot 3, 
and finally to the right lead. If dots are not on resonance, a finite bias must compensate for 
the deviation from the resonance. If another charge configuration (A^i, N2 + 1, N3) with the 
extra electron in dot 2 is also on resonance, we have four degenerate charge configurations and 
we define this point as a quadruple point (QP). If the constituent QDs are strongly connected, 
we can not use simple classical arguments for the transport and QP is defined as a point 
in parameter space where four different charge configurations have the same probabilities. 
At QPs, the extra electron can move along two alternative paths before it escapes and 
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this leads to the Aharonov-Bohm (AB) oscillationai^ in the presence of a perpendicular 
magnetic field. With empty systems, we can use single-particle picture for the transport and 
a theoretical calculation for empty TQD using transfer matrix was recently reported.^*^ For 
the transport involving states with more than one electron, the strong Coulomb interaction 
and the correlation between different charge configurations must be taken into account. A 
rate equation for the probabilities of each quantum many-body state of the system has 
been used to describe the transport through strongly interacting quantum systems weakly 
connected to the leads.— i2^>22ii2ii2^i2Sii2Ii2^i2£ In our previous work^S, it was shown that the 
interplay between the AB oscillations and spin blockade in a TQD around a QP with an 
electron trapped in dot 2 leads to spin-selective AB oscillations. Around this QP, the singlet 
states show strong AB oscillations while AB oscillations for triplet states are suppressed. 
Combined with the Zeeman splitting, these different behaviors of singlet and triplet states 
result in repeated peaks of spin-down current at lower magnetic fields and large spin-up 
current at high magnetic fields. 

In the current work, we describe in detail how we calculate the transport through a 
strongly correlated many-body system based on rate equation approach in the sequential 
tunneling regime. Higher order processes such as cotunneling will be neglected.— i^i^ We 
will focus on a TQD system which contains one or two holes. A hole is defined as absence 
of electron with respect to the fully occupied TQD with six electrons. In a resonant TQD 
where all three dots are identical, the two-hole system has spin-triplet ground state while 
two-electron system has spin-singlet ground state.— »^ The singlet-triplet gap is proportional 
to the inter-dot tunneling. In other words, the inter-dot tunneling in triangular TQD favors 
spin-singlet for two-electron system and favors spin-triplet for two-hole system. Two-hole 
system in TQD has a very interesting property that the spin structure of the system can be 
tuned by gate voltages, which is in contrast to the two-electron system where the ground state 
is always spin-singlet. By lowering the level energy of dot 2 (or raising the hole level energy 
of dot 2), two holes are localized in dot 1 and 3, and the ground state changes from triplet 
to singlet.— On the other hand, when we lower the hole level of dot 2, we can permanently 
trap one hole in dot 2 with the additional hole moving around the TQD. The four charge 
configurations with a trapped hole in dot 2 are (1,1, 2), (2, 0,2), (1,2,1), and (2,1,2) in terms of 
electron occupation numbers. As is shown later, at the QP where these configurations are 
degenerate, spin singlet and triplet states of two-hole system are degenerate. This degeneracy 
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between singlet and triplet states implies that the spin state is very sensitive to the external 
environment such as nuclear spin, and it could be used for, e.g., probing nuclear spins. We 
investigate the transport around this QP in the presence of a perpendicular magnetic field. 

The plan of the paper is as follows. In Sec. [TTl we describe our model Hamiltonian for the 
TQD, the leads, and the coupling between the TQD and the leads. The current is expressed 
using electron and hole spectral functions and master equation for the probabilities of each 
state of the system will be given. In Sec. Illlt general conditions for transport through a 
quantum many-body system in the sequential tunneling regime are discussed in terms of 
transport channels, active state and trap state. In Sec. HVl using the methods developed in 
Sec. mi numerical results in a QP with a trapped hole are presented. Negative differential 
conductance and spin-selective AB oscillations are predicted and discussed in detail. A brief 
summary will be given in Sec. |Vl 

II. METHOD 
A. Model 

We have shown previously^»i^ that the electronic properties of the TQD molecule with 
few confined electrons (A^ = 1 to 6) can be understood in the frame of the Hubbard model 
with one orbital per dot. The Hamiltonian of the TQD subject to a uniform perpendicular 
magnetic field, B = i?z, is given by 

3 1 

i=l (T i^j a i i^j 

where the operators di^ {d\^) annihilate (create) an electron with spin cr = ±1/2 on orbital i 
{i = 1, 2, 3). hia- = dl^dia- and pi = fii^ are, respectively, the spin and charge density on 
orbital level i. Each dot is represented by a single orbital with energy Sio- = ei + g*fiBBa + eo, 
with g* being the effective Lande gf-factor and fiB being the Bohr magneton. Eq is an overall 
energy shift with respect to the Fermi level of the leads which can be changed by applying 
proper voltages to external gates. Notice that Eq changes the energy differences between 
states with different number of electrons, but not the energy differences between the states 
with the same number of electrons, f/j is the on-site Coulomb repulsion of orbital i and 
Vij is the direct Coulomb interaction between two electrons in orbitals i and j. The effects 
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of the perpendicular magnetic field on the hopping matrix elements are accounted for by 
the Peierls phase factors^^i^ = tij-e^'^*'^'^ . For the three quantum dots located in the 

corners of an equilateral triangle we have (pu = 023 = 4>3i = —0/3 and 0jj = — 0ij, where 
(pB = BA/(f)Q is the number of magnetic flux quanta threading the system, with A being the 
area of the triangle and 0o = he/ e being the magnetic flux quantum. The energy spectrum 
of the isolated TQD in the presence of the magnetic field is obtained by the configuration- 
interaction method.— 

The TQD system (D) is connected to leads {r=L, R) on both sides with left lead con- 
nected to dot 1 and right lead to dot 3 (Fig. [1]). We will use non-interacting one-dimensional 
chains as our model for leads (see Appendix |X] for explicit expressions for the eigenstates 
and eigenvalues of the leads), which are connected to the TQD by 

HrD = {irD{k)cl^di^„ + h.c.^ , (2) 

(T k 

where k is the eigenstates of the lead r, iroik) is the coupling strength between the lead 
state k and dot io = 1 {io = 3) for r = L {r = R), given in Eq. (lA6p . We will assume 
that tunnel coupling tro between the lead site and the QD is very small that we can use 
sequential tunneling picture described in the following. 

B. Transport theory and master equation 

We obtain the many-body eigenstates of the TQD molecule by solving Eq. exactly for 
A^=0 to 6. The system of isolated TQD molecule and two leads is our unperturbed system, 
and Eq. ([2]) is used as the perturbation causing the transition between different states. The 
tunneling between the leads and the quantum system is assumed to be small so that we can 
consider only sequential tunneling transport. 

The net current of electrons with spin a from the lead r to TQD is 

Ir^n = i-e) i^'^^D - ^D^r) , (3) 

where ^r~>j:, is the rate of an electron with spin a moving from lead r to the TQD, which is 
given by 

5 

^r^D = E E E (^-iV.l - ^Sv) Pc^.^r («^, «7V+l) , (4) 
Af=0 ojv ojv+i 
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and ^D^r is the rate of an electron with spin cr moving from TQD to lead r which is given 
by 

5 

^D^r = E E E - (^Sv+i - ^Sv)] P-M+^^r i^N, • (5) 

Af=0 ajv ajv+l 

Here fr{£) = l/{exp[(£ — /i,,)/A;BT] + l} is the Fermi-Dirac distribution function with respect 
to the chemical potential fir of lead r and the temperature T, and Pa^ is the probability 
that the TQD is in A^-electron state with energy Ej^^. {aiy,aN+i) is the transition 
rate of the TQD state from to ajy+i by adding an electron with spin a from the lead r 
to the TQD orbital io which, by Fermi's Golden Rule, is 

k 

The transition rate from a^+i to is also given by the same {ajy, a^+i) and it describes 
the coupling between and ajy+i by the connection to the lead r. Individual transitions are 
assumed to be independent and correlation between different many body states is neglected 
in this formalism.— The current of electrons with spin a from the lead r to TQD can be 
expressed as 

5 

h' 



Af=0 k 

X [fr (4) A, (iV; zo, a- e^) - [1 - (4)] A {N + 1; ^o, a; 4) } , (7) 

where 

A,{N-t,a-e) = 27r J] P., | («^+i|4|a^) |' 5 (eJ^^^ - - e) , (8) 

We refer to Ae and Ah, respectively, as the electron and hole spectral function of the TQD 
molecule. These spectral functions contain information on the intrinsic properties of the 
TQD. In transport, what is measured is essentially these spectral functions weighted with 
the distribution functions of the leads fr and the lead-dot coupling iroik). Note that the 
probabilities of each state are determined by the master equation described below, and 
therefore the spectral functions depend on the bias. The first term in Eq. ([71) describes an 



6 



electron moving from the lead r to the TQD system increasing number of electrons in the 
TQD from to + 1. The second term describes the removal of an electron from the TQD 
with + 1 electrons back to the same lead r. 

At equilibrium, the probabilities Paj^^s are given by 

exp (-^^eZ!^) 
= ^ z'^"" ^ , (10) 

where Z is the grand partition function of the TQD molecule and /xq is the chemical potential 
at equilibrium. When a finite bias is applied between the left and right leads, the time 
evolution of the probabilities Paj^ is given by the following set of master equations. 



dP 



dt 

ajv+i r=L,R 



J2 E ^--+1 - - ^Jv)] r.(«^, «iV+l) 

r+i r=L,R 

CK]V+i r=L,R 

+ E E (^Sv - Et-.) r.(«^-i,M 

ojv-i r=L,R 

- E E [l - (i^Sv - ^Sv-.)] rr(«7V-l,a7v) , (11) 

ajv-i r=L,R 

where = '^^-^r- "^^^ ^^^^ ^^"^ third terms are the contribution from the transitions 
ctjvii — > ttAT, and the second and fouth terms are the contribution from the transitions 
aN — ^ ftAfii- For A^ = 0(A^ = 6), the summations over aAr_i(Q;7v+i) is absent in the master 
equation. If we define a vector P whose elements are Paj^ of all possible A^ and a^, the 
master equation can be cast into a matrix form: 

f = Mp.p, (12) 

where the matrix Mp is defined by Eq. (ITT]) . The initial probabilities at t = are given 
by equilibrium values Pa^^iO) = Pa^- The normalization condition that the sum of all the 
probabilities is unity, 

EE^"- = i' (13) 

N ajv 

is conserved by the master equation. 

At steady state, all the time derivatives are zero and the master equation Eq. ( fT2l) becomes 
a homogeneous system of linear equations, 

Mp ■ P = . (14) 



We can solve this equation combined with the normahzation condition Eq. (fT3l) . This 
equation has nontrivial solutions and sometimes the solution can not be uniquely determined 
with only the normalization condition. To uniquely determine the steady state probabilities, 
we need to make use of the initial conditions. Using the singular value decomposition of the 
matrix Mp, we can uniquely determine the steady state probabilities (see Appendix [B] for 
details). 

III. CONDITIONS FOR TRANSPORT THROUGH A QUANTUM MANY-BODY 
SYSTEM 

Before we move on to the numerical results, we discuss the conditions required for trans- 
port through a general quantum system connected by two leads at both ends. We will call 
the general quantum system "dot" in this section. The transport in general involves mul- 
tiple many-body quantum states and can not be reduced to a single electron description. 
Transport channels for this many-body, many-channel transport in the sequential tunneling 
regime will be defined and general conditions for transport will be given. We will assume 
zero temperature to simplify the discussion. 

Without source-drain bias (V"sd=0), the dot will be in equilibrium with the leads with 
chemical potential /xq- The number of electrons in the dot at equilibrium will be determined 
by 

fiD{N - 1) < fio < fiD{N) , (15) 
where /id is the chemical potential of the dot defined as 

fiD{N)=E^as -E% , (16) 

for the A^-electron ground state and (A^ + l)-electron ground state a^+i. The dot will 
be in the ground state of A^-electron system, a'fp . Once we apply a bias, transition from 
the ground state to other states will occur. By succession of such transitions, electrons 
can move from one lead to the other lead. During this process many of the dot states will 
be accessed and on the average we can assign a probability that the dot will be in each 
many-body state. Probabilities will change in time initially, but reach stationary values at 
steady state. 
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To be specific, let us assume that we apply a bias Vsd between the two leads such that 
A^L = yUo + and yUj? = yUo — e\4d/2. For a transition from an A^-electron state to an 

(A^ + l)-electron state Pn+i by adding an electron in the dot from lead r = L, R to occur, 
the incoming electron must have the energy of E^^^^ — E^^. Electrons in the lead r have 
energies below the chemical potential /i^ and therefore the condition E^^^^ — E^^ < /i^ must 
be satisfied. In addition, the states aN and I3n+i must be connected by adding an electron 
from the lead r, i.e., r^(a7v, /^at+i) 7^ 0. Similarly, a transition from an (A^ + l)-electron 
state /^AT+i to an A^-electron state a at by moving an electron from the dot to the lead r is 
allowed if E^^^^ — E^^ > /i^ and Tr{aN, Pn+i) 7^ 0. In this case, there must be an empty 
state with energy E? — E^ in the lead r. We define a set of two successive transitions 

as a transport channel if [i) these transitions are allowed and [ii) the initial state or 
Pn+1 is either the ground state or accessible from the ground state by successive allowed 
transitions. The L or i? on top of the arrows means the transition is allowed by connection 
to the left or right lead. If two states are involved in a transport channel such as — >■ 
Pn+1 C(N, we will simply call the pair {un, Pn+i) a transport channel. We will define any 
state that participates in one or more transport channels an active state. In most cases, if 
there are one or more transport channels, the current flows through the system. Exception is 
when there is a trap state. A state •yjy is defined to be a trap state if (i) •jn is accessible from 
the ground state through successive allowed transitions and (ii) the transition 77V aAr±i is 
not allowed for any active state aN±i. When we apply a bias to the system, the probability 
of the trap state monotonically increases in time and in the steady state condition, the 
system will be completely in the trap state and current can no longer proceed. The (1,1) 
triplet state in DQD spin-blockade system^^ with forward bias is an example of trap state. 
To summarize, the general conditions for transport via sequential tunneling are (i) there is 
at least one transport channel and {ii) there is no trap state. 

We can understand the expression for total current in Eq. ([7]) in terms of transport chan- 
nels, active states and trap states. If there is no trap state, then each active state will have 
a finite probability. The electron spectral function [Eq. ([8])] is sum over all transitions from 
an active state aN to any state oat+i. For the total current, this electron spectral function is 
weighted by the Fermi-Dirac function. At zero temperature, this weighted spectral function 
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is sum over all transitions only between active states. Similarly, the hole spectral function 
is sum over all transitions from an active state ajv+i to any A'"-electron state. It is weighted 
with I — fr and the weighted hole spectral function is sum over all transitions between active 

states. If there is a trap state 77V, then only P^j^(=l) is nonzero and the weighted spectral 
function is sum over transitions from 77V to other active states. Since there is no allowed 
transition from to any active state by definition of the trap state, the weighted spectral 
function is zero and there is no current. 

IV. RESULTS 

We use the effective Rydberg Ry and the effective Bohr radius as of the host semicon- 
ductor material as our units for energy and length, respectively. For GaAs, Ry = 5.93 
meV and = 9.79 nm. The three QDs are located at the vertices of a equilateral trian- 
gle and the inter-dot distance is 6.25. The magnetic field is measured with the number of 
flux quanta 0b through the triangle. We use parameters Ui — U — 2.5, Vij — V — 0.5, 
tij — —t — —0.05, for all i and j. The tunnehng strengths ti and tu in the leads are 
assumed to be large so that the density of states of the leads are nonzero in a wide range 
of energy. We use ti = = —4.0 in the calculation. The current is calculated in units 
of Iq = eltLDl"^ /h\tL\, and conductance in units of Gq = e^|tLDp/^|ti|. We assume that 
thD — tRD and that they are small enough to justify the sequential tunneling picture. In 
these effective units, the current and the conductance do not depend on the dot-lead tun- 
nehng tiD and tuD- The effective g'-factor is taken to be —0.44. The chemical potentials of 
each leads are fiL = A*o + and fJ>R = fJ^o — eV^d/2. We set the equilibrium chemical 

potential /xq to be zero. All the calculations are done at temperature T=0.001. 

A. Quadruple point with a trapped hole 

We consider a quadruple point where a hole is trapped in dot 2. A hole in the TQD is de- 
fined as the absence of an electron with respect to the fully occupied state with six electrons. 
At this quadruple point, the four resonant charge configurations are (2, 1,2), (1,1, 2), (2, 0,2) 
and (2,1,1). If we neglect the tunneling between the constituent dots of the TQD molecule. 
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the QP can be determined by equating the energies of the four charge configurations, 



E(2, 1, 2) = E{1, 1, 2) = Ei2, 0, 2) = E(2, 1, 1) = Eqp . (17) 

With El = = e and E2 = e + A, the energies of the four configurations are 

E{2, 1,2) = 5E + A + 2U + 8V, (18) 
E{1, 1, 2) = E(2, 1, 1) = 4e + A + f/ + 51/ , (19) 
E{2, 0, 2) = AE + 2U + AV , (20) 

and from Eq. f|T7|) we obtain e = — f/ — 31^ and A = f/ — At this classical QP the level 
energies of the QDs are 

E, = E3 = -U-3V , (21) 

82 = -W , (22) 

and the total energy is 

Eqp = -2U - 8V . (23) 

When we have tunneling between the QDs, the charge configurations in terms of the pop- 
ulation of each individual dots are not eigenstates of the coupled TQD molecule system. 
The quadruple points are defined then as the point where all the resonant charge config- 
urations have the same probabilities, which can be determined numerically. The QP with 
a trapped hole where the four charge configurations (2, 1,2), (1,1, 2), (2, 0,2), and (2,1,1) have 
same probabilities corresponds to the level energies Ei = E3 = —4.05 and E2 = —2.095 with 
£0 = 0. Notice that these values are quite close to the values of the classical QP because the 
tunneling t is relatively small compared to the Coulomb interaction parameters U and V. 

In terms of hole occupation numbers, the charge configurations (0,l,0)/i, (l,l,0)/i, (0,2,0)/^, 
and (0,1,1),, are degenerate where (iVi, N2, N3)h={2 - Ni,2 - N2, 2- N3). The lower part 
of the energy spectrum of the TQD at this QP is shown in Fig. [2]^a) . The two thin black 
dotted lines are the energy levels of spin-up and spin-down states of the trapped hole. There 
is no AB oscillations for the single-hole states because it is localized in dot 2. For two- 
hole states, neglecting very high energy states consisting mostly of configurations (2,0,0),,, 
(0,0,2)/j and (1,0,1)/,, we have nine two-hole states (3 singlet and 6 triplet states). Two- 
hole singlet states (red solid curves) show strong AB oscillations since they consist of three 
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configurations (1,1, 0)^, (0,2,0)/i, and (0,l,l)/i. That is, there is one hole trapped in dot 2 and 
an extra hole moves around the three dots. Two-hole triplet states mainly consist of only 
two configurations (l,l,0)/i and (0,1, 1)/^ because (0,2,0)/j configuration can be only singlet 
due to the exclusion principle. Therefore, triplet states (blue dashed curves) show only small 
AB oscillations which can be ascribed to higher energy triplet configuration {1,0, l)h. This 
suppression of AB oscillation in triplet states also occurs at the QP with a trapped electron 
in dot 2.— The energy spectrum of the TQD with a trapped hole has two main differences 
compared to the TQD with a trapped electron. First, at zero magnetic field, two- hole spin- 
singlet and triplet states are almost degenerate. Second, the phase of the AB oscillation of 
the singlet state is shifted by vr. Without Zeeman splitting, these differences would lead to 
7T phase-shifted AB oscillations with less amplitudes in current because the triplet transport 
channel will have finite contribution with very small AB oscillation. With Zeeman splitting, 
the single-hole and two-hole ground states are on resonance only at zero magnetic field. We 
can tune the relative energy difference between two-hole {i.e., four-electron) and single-hole 
{i.e., five-electron) states by changing the overall energy shift Eq. The two-hole states obtain 
Aeq while single-hole states obtain 5£:o- Therefore, single-hole states gain additional energy 
6q with respect to two-hole states. It is easier to understand how the transport channels 
contribute to the total current if the one-hole and two-hole states are energetically separated. 
Figure[2t^b) is the energy spectrum with eo=-0.15, where the single-hole states are well below 
the two-hole states. Interesting phenomena such as negative differential conductance and 
spin-selective AB oscillations occur in this system. We will present results for this system 
at various biases and magnetic fields. 

B. Spectral functions and transport channels 

To explain how the current is related to the spectral functions and the transport channels, 
we choose a specific case with 0^=3.3 and eV^d=0.3. The system is schematically given in 
Fig. [31 The single-hole states with spin up (/5|) and down {P^) are the lowest levels. For 
two-hole states, as is spin singlet and is spin triplet with total ^^=±1, respectively. 

Solid black (red) arrows represent transitions from two-hole states to single-hole states by a 
spin-down (spin-up) electron moving from the left lead to the TQD system. These transitions 
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satisfy 

E^-E^<fiL = f^o + ^ = 0.15 , (24) 

and rL{<y,(3) 7^ 0. ^ and other higher energy levels are outside the transport window 
and are not accessed during transport. Notice that the transition a!^^ —>■ jS^ is energetically 
possible but not allowed by spin-blockade since adding a single electron cannot change the 
total spin 5*^ from +1 to -1/2. The transition from to /5| can also occur by adding an 
electron from the right lead since the incoming electron must have energy -0.17, which is 
below the chemical potential of the right lead /i/j=-0.15. This is the only transition that 
transports an electron from the right lead to the TQD, and is represented as a dotted red 
arrow pointing downward. Dotted black (red) arrows pointing upward represent transitions 
from single-hole states to two-hole states by a spin-down (spin-up) electron moving from the 
TQD system to the right lead. These transitions satisfy 

-E^>f,n = l^o-^ = -0.15 , (25) 

and Tji{a,(3) 7^ 0. There are four pairwise transport channels (05, /5^), (0:5, /5j^), (a^\/?j^), 
{a!^\(3i) as well as other transport channels such as — >• 05 — etc. There is no trap 
state in this case. 

The contribution of each transition to the current can be understood by spectral functions. 
Figure H] shows the spectral functions for this system. Due to the symmetry of the system 
and the assumption tLD=^/?D, spectral functions for i=l and i = 3 are the same. Spectral 
functions are sum of delta-functions and the coefficients of each delta-function are plotted 
here, (a) is the electron spectral function and (b) is the hole spectral function. Electron 
(Hole) spectral function shows transitions from any active two-hole (single-hole) state to 
any single-hole (two-hole) state, not necessarily only to active states. But not all these 
transitions actually occur during the transport. Only those transitions within the transport 
window occur and contribute to the current. For the current from the left lead to the TQD, 
the summation is over the electron spectral function weighted by the Fermi function of 
the left lead. Black (Red) columns in (c) show the weighted electron spectral function for 
adding a spin-down (spin-up) electron, which correspond to the solid black (red) arrows in 
Fig. [31 Transitions aip^ (3^ and a^^ (3i are induced by spin-down electrons with same 
energy and the highest black column represents sum of these two. This weighted electron 
spectral function is the same with the bare spectral function in the energy range shown here 
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because all the transitions from active two-hole states to any of the and are in the 
transport window. For the current from the TQD to the left lead, the hole spectral function 
is weighted by I-Jl [(d)]- This weighted hole spectral function is zero since no transition 
from a single-hole state to a two-hole state emits an electron with energy larger than the 
chemical potential of the left lead (there is no solid arrows pointing upward in Fig. [3]). For 
the current from the TQD to the right lead, we need to sum over the hole spectral function 
weighted by 1 — Black (Red) columns in (f ) correspond to the dotted black (red) arrows 
pointing upward in Fig. [3l For the current from the right lead to the TQD, the electron 
spectral function is weighted by //? [(e)], which has only one small peak corresponding to 
the transition a^^^ to /9| represented by the dotted red arrow pointing downward in Fig. O 
These weighted spectral functions show the transitions between active states {i.e., transitions 
shown in Fig. [3]) and only these transitions form transport channels and contribute to the 
current. The net current from the left lead to the TQD Il^d can be obtained by summing 
the peaks in (c) and the net current from the TQD to the right lead Id-^r is the sum of the 
peaks in (f) minus the peak in (e). The net currents Il^d and Id^r are the same, satisfying 
the steady state condition. 

C. Negative differential conductance 

When we increase the bias, the transport window expands and more transport channels 
are involved in the transport. Whenever a new transport channel is introduced, the prob- 
abilities of each states must be redistributed and the spectral functions and the current 
change accordingly. In the range of biases where no new transport channel is introduced 
the current remains fiat and therefore nonzero differential conductance signifies introduc- 
tion of new transport channel. Figure shows how the currents[(a) and (d)], differential 
conductances [(b) and (e)], and the probabilities[(c) and (f)] of each states change for the 
same system as in Fig. [2](b) as we increase the bias at two different magnetic fields 0_b=O 
and 3.3. At zero magnetic field [(a)~(c)], the single-hole state is doubly spin-degenerate and 
the current changes when a two-hole state enters the transport window as we increase the 
bias. At each current plateaus the probabilities of each active states are the same, which 
results from the assumption that tLZ)=^i?D- The differential conductance peaks at the energy 
differences between two-hole states and the single-hole ground state. At finite magnetic field 
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0B=3.3 [(d)~(f)], the spin-degeneracy of the single-hole states is lifted and we have multiple 
A^- and (A^ + l)-particle states. The current changes whenever additional transport channel 
is introduced by increasing bias, which is more complicated than at zero magnetic field since 
the two single-hole levels have different energies. The probabilities of each active states at 
current plateaus are not the same in this case. In most cases the additional transport chan- 
nel results in the increase in current, but sometimes it leads to the decrease in the current as 
is the case for the decrease in current at eVsd — 0.7 for (f)B = [Fig. [5](a)], and at eVsd — 0.62 
and 0.69 for 0^ = 3.3 [Fig. EJ^d)]. This negative differential conductance can be explained 
as a result of the interplay between the different coupling strength r.r(a,/?) for singlet and 
triplet states and the probability redistribution of many-electron states with increased bias. 
The single-hole states are more strongly coupled to the triplet two-hole states than to the 
singlet two-hole states (F^ is larger for the triplet) because the triplet states consist of two 
configurations [(l,l,0)/i and (0,1, 1)/^] while the singlet states consist of three configurations 
[(l,l,0)/i, (0,2,0)/i, and (0,l,l)/i], and the configuration (0,2,0)^ is not connected to the single 
hole configuration (0,l,0)/i by adding or subtracting a hole from the leads. If increasing the 
bias makes the singlet state above the triplet states become an active state, we have more 
transport channels. But the probabilities of all the active states will be redistributed and the 
triplet active states will have less probabilities. If the current decrease due to the reduction 
of the triplet probabilities is larger than the current increase due to the introduction of new 
transport channels involving the singlet state, the total current decreases as we increase the 
bias. 

Figure M shows the electron spectral functions weighted by /l before [(a) and (c)] and 
after [(b) and (d)] the current decrease. The electron spectral functions weighted by l-fi 
are all zero for this case, and the current is simply proportional to the sum of the heights 
of all columns in the figure. At zero magnetic field [(a) and (b)] we show only spin-down 
component because spin-up component is the same due to the spin degeneracy. Before the 
current decrease, at eV"sd=0.6 [(a)], the two high columns correspond to the transitions from 
spin-triplet states to single-hole states and the two small columns which almost overlap with 
one of high column are from the transitions from spin-singlet states. The three columns are 
very close together due to the singlet-triplet degeneracy of the two-hole system. After the 
current decrease, at eV^d=0.8 [(b)], there is additional peak at around e ~ —0.35, which 
correspond to the new allowed transitions involving the higher singlet state. Notice that the 
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decrease of the two peaks of triplet transitions is bigger than the new singlet peak, which 
leads to the negative differential conductance. At finite magnetic field [(c) and (d)], spin 
degeneracy is lifted and there are more peaks in the weighted spectral function. The higher 
peaks correspond to transitions from triplet states and lower peaks correspond to transitions 
from singlet states. At eV^d=0.8 [(d)], the two new peaks correspond to the transitions from 
the third singlet state to the two single-hole states. Once again, the peaks of triplet transport 
channels decrease and the total current decreases. 

D. Spin-selective Aharonov-Bohm oscillations 

The different oscillatory behavior of singlet and triplet states in the magnetic field can lead 
to spin-dependent transport phenomena. Left panels of Fig. [7] [(a)~(c)] show the current, 
current polarization, and differential conductance as a function of the magnetic field at 
eKd=0.25 for the same system as in Fig. [2](b). At this bias, the lowest singlet two-hole state 
(as) can form transport channels with both the spin-up (/5|) and spin-down (jSi) single-hole 
state at lower magnetic fields. (3^ and have the same probabilities in this regime and 
the spin-up and spin-down components of current are the same and the spin-polarization 
of the current is zero. The AB oscillation makes the singlet state as oscillate in and out 
of the transport window, which leads to the oscillations in current and conductance. At 
higher magnetic fields, the spin-triplet state with Sz=+l (a^"*) comes into the transport 
window with respect to and Since /3| has S'^=-l/2, it cannot form a transport channel 
with due to spin blockade. Therefore, the transport channel {a^\ jS-^) which carries 
spin-down current dominates the transport at high magnetic fields. At higher bias, with 
eVsd=0.3 [Fig. [7] (d)~(f)], the pair of states {a^\ (3^) always forms a transport channel and 
spin-down current is dominant at all magnetic fields. The singlet state as can form transport 
channels with both /?| and Pi states and contributes equally to the spin-up and spin-down 
currents. The formation of transport channels involving spin-singlet as leads to probability 
redistribution and decreases the probability of the triplet state a!^\ Therefore, we have 
higher current polarization when the singlet state as does not participate in the transport 
and vice versa. 

Figure [S] summarizes both the negative differential conductance and the spin-selective 
AB oscillations at finite bias for the system in Fig. Mjo). (a) and (b) show the spin-down 
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and spin-up component of the differential conductance as functions of magnetic field and the 
bias, which clearly show different AB oscillations for spin-down and spin-up conductance, 
(c) is the total differential conductance and the dark trace at large bias due to introduction 
of transport channels involving the high energy singlet state indicates the current decrease 
for increasing bias. Since nonzero differential conductance indicates introduction of new 
transport channels and the bias for new transport channel corresponds to the energy dif- 
ference between the two constituent states of the new transport channel, the differential 
conductance can be used as a spectroscopic tool. Comparing Figure [8](c) with the energy 
spectrum Fig. [2](b), we can see the resemblance. Differences are that singlet levels are split 
into two and three Zeeman-split triplet levels give only two nonzero traces in differential 
conductance G. The lowest singlet state forms transport channels with (3^ and /3| at the 
same bias, as we mentioned earlier, and therefore leads to only one nonzero trace in G, while 
higher singlet states can form transport channels with (3^ and (i^ at different bias and hence 
two nonzero traces in G. For the triplet states, triplet states with S'^=+l (-1) can not form 
transport channel with (3i {(3^) due to spin blockade and therefore we have only two nonzero 
traces in G for three triplet states. 

To compare with the spin-selective AB oscillations of TQD system with a trapped electron 
in the linear response regime,— we consider the TQD slightly off the exact QP with a trapped 
hole. Figure [9] shows the spin- selective AB oscillations with Eq = —0.025 in the linear 
response regime where eV^d = 1-0 x 10~^. The spin-up single-hole state crosses the singlet 
two-hole state several times until it crosses with the triplet state at higher magnetic field. 
The repeating conductance peaks are for spin-up current and the large peak at strong enough 
magnetic field is for spin-down current. Compared to the case at the QP with a trapped 
electron, the behaviors of spin-up and spin-down currents are reversed since the transport of 
a hole with spin-up (spin-down) corresponds to the transport of an electron with spin-down 
(spin-up) in the opposite direction. The phase of the oscillation is also shifted by vr due to 
the TT phase shift of oscillations of the singlet states. The current is ultimately suppressed 
at very high field where the triplet state is the ground state and spin-up single-hole state is 
outside of the transport window. 
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V. CONCLUSIONS 

We presented a theory of the tunnehng transport through a TQD around a QP with a 
trapped hole in dot 2 where spin-selective AB oscillations occur. A detailed description of the 
formalism for the transport calculation and general conditions for transport through multiple 
many-body states were given. It was shown that the interplay between the introduction of 
new transport channels and the probability redistribution can lead to negative differential 
conductance and that the differential conductance can be used as a spectroscopic tool. The 
transport in the magnetic field is sensitive to the spin of the carriers and the spin structure 
of the TQD system in a triangular geometry due to the strong Coulomb interaction and the 
interference effects. This spin-selective AB oscillations show different behavior in the linear 
response regime and with a finite bias. The TQD system with singly connected leads can be 
considered as a minimal quantum dot network (QDN) circuit and the formalism developed 
here can be used for general QDN circuits. The multi-channel transport can be analyzed 
using the spectral functions of the QDN and transport channels formed by many-body states. 
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APPENDIX A: NONINTERACTING ONE-DIMENSIONAL CHAIN 

The leads are described by the Hamiltonian of non-interacting one-dimensional chains 



where Sq is the level energy of each site and tr is the tunneling between sites in lead r = L. R. 
m is from — A''^ to —1 for the left lead and from 1 to A^ for the right lead. The eigenstates 
of this tight-binding chain is 





(Al) 
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where k = 2mi/ (Naa) with n = 0,1,2, Na — 1. The eigenvalues are spin degenerate and 
given by 

el = el + 2tr cos ka , (A2) 
and the density of states per site for each spin is 

PrAe) = . |, .\ , e i2\tr\ -\e-el\) , (A3) 

where 9 is the step function and k^ is determined by e = eg + 2t,. cos(A;ea). 
The TQD and the leads are connected by tunneling Hamiltonian 

HrD = {trDci^^Ao^ + h.c) , (A4) 

a 

where mo (—1 for r = L and 1 for r = R) and iq (1 for r = L and 3 for r = R) are the two 
adjacent sites of the lead and the dot connected by the tunneling and Ud is the tunneling 
element connecting the two sites. Using the eigenstates of the lead chains, we obtain 

^^•^ = {^'■^^^''^kadioa + h.c^j , (A5) 

cr k 

where 

imik) ^ . (A6) 

APPENDIX B: SOLUTION OF THE RATE EQUATION AT THE STEADY 
STATE 

In this appendix, we present a method to find the steady state solution of the rate 
equation. The master equation in matrix form is give by 

^ = Mp.p, (Bl) 

where the matrix Mp is defined by Eq. (11 II) . The initial probabilities at t = are given by 
equilibrium values 

exp f_<vz^) 
Pa AO) = PZ = ^ /"^ ^ , (B2) 

For stationary cases {t — ^ oo), dPaj^/dt = for all a^. Then we obtain a system of linear 
equations for a^r 

Mp ■ P = . (B3) 
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This equation by itself does not uniquely determine Pq,^ because the matrix Mp is singular. 
We can see this by the fact that the summation of all the elements of Mp is zero and therefore 
the set of equations in Eq. (1B3P are not linearly independent. We need more conditions to 
uniquely determine the steady state solution. One constraint is the normalization condition 
of the probabilities: 



If all the states take part in the transport, then Eqs. (]B3|) and (]B4|) would uniquely determine 
the steady state probabilities. But if some states don't participate in the transport, we need 
more conditions (If we consider Eq. (1B3P as an eigenvalue problem with eigenvalue 0, the 
eigenvalue can be degenerate). Considering that the master equation gives unique solution 
with the initial condition, we need to make use of this initial condition to solve steady state 
solution in this case. 

For this, let us use the singular value decomposition (SVD) of Mp. 



where U and V are orthogonal matrices and D is a diagonal matrix. The diagonal elements 
of D are called singular values. All singular values are zero or positive and we assume that 
the diagonal elements of D are in descending order. The number of zero singular values 
(A^o) is the dimension of the null space of Mp, and the last Nq column vectors of V form 
a basis set for the null space. If only one singular value is zero, the last column vector of 
V defines the one-dimensional null space of Mp and the the normalization condition would 
be enough to uniquely determine the probabilities. But, in general, there can be more than 
one singular values. For Np x Np matrix Mp, Np = A^p + A"o where A^p is the dimension 
of the range and Aq is the dimension of the null space. Eq. ( IBll) can be written as 




(B4) 



Mp = UDV 



■T 



(B5) 



dz 
'dt 



Dy , 



(B6) 



where 



z = 



y = 



v^p . 



(B7) 
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Separating the null space (B) and its complementary space (A), it becomes 



d 
di 




is a diagonal matrix whose diagonal elements are the nonzero singular values of Mp. 
Since time derivative of zp is zero, we obtain 

ZB(t) = zb(0) (B9) 

at all time t. At t —>■ oo, all time derivatives are zero and we obtain 

^J^(LpSl = D.y.(« - oo) = (BIO) 
^ y^(i^oo) = (BU) 

since is a diagonal matrix with all nonzero diagonal elements. Using 

y = V^P = V^Uz = Wz , (B12) 

^-^-U-) , (B13) 
^Wba Wbb J \zb J 

we get 

yA(oo) = WaaZa(oo) + WabZb(oo) = (B14) 
and we obtain an equation for Zyi(oo) 

WaaZaM = -WabZbM (B15) 
^ z^(oo) = - (WAAr' W^bZb(O) , (B16) 

where we used the result that z^ is constant in time. We can prove that (W^^)~^ exists 
as follows. Since Mp has a null space of dimension A^^q, Mp has eigenvalue of degeneracy 
A^o- From Eq. ( IBSll . we have 

U^MpU = DW = I I . (B17) 
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Since U'^MpU has also eigenvalue of degeneracy A^o, 

det (U^MpU - AI) = A^»/(A) , (B18) 

where /(A) is a polynomial with /(O) 7^ 0. Using Eq. (]B17p . 

det(DW-AI) = A^Met (DaW^a - AIaa) (B19) 

and A = is not a solution of det (D^Waa — AIaa) = 0. Therefore, 

det (DaWaa) = detBA detW^A 7^ , (B20) 

and, since is a diagonal matrix with nonzero elements, we obtain detW^A 7^ and 
(Waa)~^ exists. Once we have z(oo), we can find P at steady state using 

P(t 00) = Uz(oo) . (B21) 

The normalization condition is automatically satisfied in this method, because the master 
equation conserves the normalization. 
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FIG. 1: Schematic diagram of a triple quantum dot molecule connected to leads, fs are the 
tunneling matrix elements that connect two different sites. Each lead is modeled by a semi-infinite 
noninteracting chain which is connected to the TQD at one end. 
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1 2 3 4 5 6 



FIG. 2: (Color online) (a) Energy spectrum at the QP where a hole is trapped in dot 2. Four charge 
configurations (0,l,0)/j, {l,l,0)h, (0,2,0)/i, and (0,l,l)/i in terms of the hole occupation numbers have 
the same probability at zero magnetic field. Black dotted curves are for single-hole states with spin 
up and down. Red solid (Blue dashed) curves are spin-singlet (spin-triplet) states for two-hole 
systems, (b) Energy spectrum with the overall shift eo=-0.15. By changing eq, the single-hole 
states move upward for positive Eq and downward for negative Eq with respect to the two-hole 
states. 
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FIG. 3: (Color online) Schematic picture describing the transitions at (/)b=3.3 and eT4d=0-3- /3's 
are single-hole levels with spin up and down, and a's are two-hole levels, as is spin singlet and 
a^''^'' are spin triplet. Solid (Dotted) arrows are transitions allowed by connection to the left 
(right) lead. Transitions represented by downward arrows are from a two-hole state to a single- 
hole state by adding an electron from the lead, and upward transitions are from a single-hole state 
to a two-hole state and emit an electron to the lead. The black (red) color represents that the spin 
of the electron transporting in the transition is down (up). Energies and the probabilities of each 
level are also given for reference. 
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FIG. 4: (Color online) Spectral functions at (/)b=3.3 and eT4d=0.3. Black (red) columns are 
for a=l (cr=|). (a) Electron spectral function, Ae{N=A;i,a;e). (b) Hole spectral function, 
Ah{N=5;i,a;e). i=l or 3 gives the same spectral functions since we have symmetric system 
and assumed tLD=tRD- (c) and (e) are electron spectral functions weighted by the Fermi function 
fi and fa respectively, and represent allowed transitions from a two-hole state to a single-hole state 
by an electron moving from the lead to the TQD system, (d) and (f) are hole spectral functions 
weighted by I-Jl and respectively, and represent allowed transitions from a single-hole state 
to a two-hole state by an electron moving from the TQD system to the lead. No transition moves 
an electron from the TQD to the left lead and therefore there is no visible peaks in (d). These 
weighted spectral functions show all the transitions during the transport depicted by arrows in 
Fig. El 
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FIG. 5: (Color online) Currents, differential conductance, and probabilities as functions of the 
bias at two different magnetic fields 0b=O [(a)~(c)] and (/)b=3.3 [(d)~(f)], for the same system 
as in Fig. EJb). For the currents and conductances, black solid curves are the total currents or 
conductances and red dashed (blue dotted) curves are for spin down (up) components. For the 
probabilities, black dotted curves are for single- hole states and red solid (blue dashed) curves are 
for spin-singlet (spin-triplet) states. 
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FIG. 6: (Color online) Electron spectral functions weighted by fi before [(a) and (c)] and after [(b) 
and (d)] current decrease. At zero field [(a) and (b)], spin-down and spin-up components are the 
same due to the spin degeneracy. At finite field [(c) and (d)], more peaks appear due to the lifting of 
the spin degeneracy. The introduction of the transport channels involving the high-energy singlet 
state as the bias increases leads to redistribution of the probabilities and reduction of the current 
through the triplet transport channels which are dominant channels. The total current decreases 
since the decrease through the triplet transport channels is bigger than the increase through the 
new singlet transport channels. 



30 



eV =0.25 

sd 



eV =0.3 

sd 




FIG. 7: (Color online) Spin-selective AB oscillations at two different values of bias, eT4d=0.25 
[(a)~(c)] and 0.3 [(d)~(f)]. Top panels show the spin-down (black solid curves) and spin-up (red 
dashed curves) currents. Bottom panels show the differential conductance of each spin species. 
Middle panels show the spin polarization of the current, defined as pol = (I^ — /|)/(/| -|- /|). 
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FIG. 8: Differential conductance as a function of the magnetic field and the bias for the system in 
Fig. [2|^b). (a) and (b) show spin-down and spin-up components separately and (c) shows the total 
differential conductance. 
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FIG. 9: (Color online) Spin-selective AB oscillation around a QP with a trapped hole in dot 2 
(eo = —0.025). Spin-up conductance Gu shows repeated peaks at lower magnetic fields and spin- 
down conductance Gd has large peaks at large magnetic fields. Inset shows the energy spectrum 
which is shifted from the QP by eo = —0.025. 
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